Phases of a fermionic model with chiral condensates and Cooper pairs 

in 1+1 dimensions 
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We study the phase structure of a 4-fermi model with three bare coupling constants, which po- 
tentially has three types of bound states. This model is a generalization of the model discussed 
previously by A. Chodos et al. [Phys. Rev. D 61, 045011 (2000)], which contained both chiral 
condensates and Cooper pairs. For this generalization we find that there are two independent renor- 
malized coupling constants which determine the phase structure at finite density and temperature. 
We find that the vacuum can be in one of three distinct phases depending on the value of these two 
renormalized coupling constants. 
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I. INTRODUCTION 

Spontaneously broken symmetry phases occur in high- 
energy and condensed matter physics by varying the tem- 
perature, density, or a parameter controlling the interac- 
tion. These include chiral and color broken symmetries 
in high-energy physics (sec e.g. Refs. 0,0,0), as well as 
superconducting, and spin and charge density waves or- 
dered phases in condensed matter systems Q . The phase 
diagram of these different physical systems depends in 
general on the compatibility of the different broken sym- 
metries and usually is studied in the framework of the 
Ginzburg-Landau- Wilson energy functional |f|- 

Recently we have studied the phase diagram of a rela- 
tivistic model with both chiral condensates and Cooper 
pairs 0, as well as a non-relativistic model with ferro- 
magnetic superconductivity 0- For ferromagnetic su- 
perconductors it is known y| that for spatially inhomo- 
geneous order parameters the ferromagnetic and super- 
conducting phases coexist in a small region of the phase 
diagram. In the context of relativistic field theory this 
has been recently studied by Rajagopal's group 0. 

In this paper we study a simple relativistic model 
in (1+1) dimensions which combines the Gross-Neveu 
model with a model for Cooper pairs P ms an 
interaction which splits the masses of the fermions. The 
symmetric version of this model has been discussed pre- 
viously [|| , and exhibits (at the mean- field level) a phase 
diagram which mimics some of the features expected for 
QCD with two light flavors of quarks. The model has 
a well defined l/N expansion and is asymptotically free 
so that it does not suffer from the cutoff dependencies 
of 3 + 1 dimensional effective field theories considered by 
others HE El. The question we are addressing here 
is whether the splitting of the fermion masses can lead 



to phase coexistence in 1+1 dimensions. What we find 
is that no such phase coexistence occurs for static mean 
fields. 

The presence of finite temperature condensates in 
one spatial dimension violates the Mermin- Wagner theo- 
rem 0| . Nevertheless, Witten has argued 01 that it is 
still meaningful to study the formation of such conden- 
sates to leading order in l/N, as the large-N expansion 
can give qualitatively good understanding of the correla- 
tion functions, even when it gives the wrong phase tran- 
sition behavior. The quality of the approximation should 
improve as we increase the number of spatial dimensions 
and the mean-field critical behavior becomes exact above 
3+1 dimensions. Therefore, the 2-dimensional realization 
of the model presented in this paper is a "toy" model ex- 
hibiting some of the features expected to be true in 3+1 
dimensions such as the restoration of the symmetry at 
high temperatures. 

The paper is organized as follows: In Sec.[n]we discuss 
the model and derive the (1+1) dimensional effective po- 
tential in the leading-order large N approximation. The 
rcnormalization of the effective potential is performed in 
Sec. IIIII We discuss the possible phases of the vacuum 
state in Sec. IIVI and study the phase structure of our 
2-dimcnsional model in Sec.0 Our conclusions are pre- 
sented in Sec. IVII 
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II. MODEL 



We would like to make the simplest generalization of 
our previous model which has N families of a two- 
flavor model, but also contains an interaction that treats 
the two flavors differently. The flavors can be thought 
of as up and down quarks, or proton and neutron. Thus 
we are lead to consider a model described by the 2N - 
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W°'' a) ] 2 



component Lagrangian 
C = ^•°)(i$-//7 )V> (< ' a) 

+ ^[^(0^(*.«)][^a)^a.«)] + 2t 

The superscripts indicate the flavor indices, i = 1,2 and 
the family index a = 1, 2, . . . , n . We will treat this model 
in the leading order in largc-N approximation which is 
a mean-field approximation. The coupling constants gf 
and G 2 must generically scale as X/N, with A fixed, in 
order to obtain the large TV limit. For simplicity, in the 
following we drop the index a. 

To perform the Hubbard-Stratonovich transforma- 
tion we first introduce auxiliary fields by adding 
to the Lagrangian the terms 

[^ (4 V 4) ]} 2 (2) 



AC = - 



9ll 



5 L{ m2 + 52 2 [^)r 3 , i ^ W ]} 2 
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In C = C + AC, the terms quartic in the fermion fields 
cancel, and we can formally write 



a = 



i 



~G r 

-h 2(e5t) 



(3) 



-2(eB) 



where e = — i(T 2 , and /i = 7°(i^ — /try — md — m 2 7°T 3 . 
We integrate out ^ and t/^ and realize that the resulting 
action is proportional to N. Performing the integration 
over the auxiliary fields by steepest descent and Legendre 
transforming the generating functional we obtain, in the 
standard manner, the leading order in large-N effective 
action 



T cS (m 1 ,m 2 ,M) 
with M 2 = AB^B~h 



(4) 
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AG 2 ) 



a 2 h<T2 



V^( mi ,m 2 ,M), 
and [H 



V<V = ^Trln{h T h) + ^Trln[l + M 2 /{hh T )] . (5) 

The gap equations arc obtained by setting to zero the 

derivatives of the effective action with respect to mf 
and M 2 , i.e. 



S m 'T = 



<5 M ^r = 



1 

2ffF 



[d d k] { 



Tr 



1 8{h T h) 



Af 2 Tr 



f 



h T h Sm 2 
1 S(hh T ) 



1 

4G 2 



Af 2 

[d d fc] Tr 



hh T hh T 5m 2 
1 



M 2 + ^ T ' 



(6) 

}■ 

(7) 



with [d d k] = d d fc/(27r) d . The effective potential is ob- 
tained as 

m 2 m 2 Af 2 

V eS (m h M)=^ + g| +— +y«(m 1 ,m 2 ,M) . (8) 

In (1+1) dimensions, a convenient representation for 
7 M is given by the Pauli matrices, as 7 = <j\ and 7 1 = 
—io 2 . This gives 

h = (k - (i) - Tfi\U\ - fci(j 3 - m 2 criT 3 , (9) 
h T = - (fc + fi) - m\0\ + fcitj 3 — m 2 criT 3 . (10) 

The pure Cooper pair version of this model, mi = m 2 = 
0, has been discussed in Ref. ^lj, while the symmetric 
version (m 2 = 0) was studied in Ref. [||. Here, we follow 
closely the approach outlined in these references. 



A. symmetric case: m,2 = 

In the symmetric case @, the effective potential is 
isospin independent, and the isospin trace results in a 
multiplicative factor of 2. 

Considering first the zero temperature case we intro- 
duce the notations: h T h = A+B-a and hh T = A'+B'-a. 
Here, we have 0: A 2 - B ■ B = A' 2 - B' ■ B'. Hence, 
we obtain the derivatives 



1 (A> + M 2 )6 mi A> 



B' -8 m 2B' 



mi (k 2 - k 2 + )(k 2 - fc 2 ) 
2(A' + M 2 ) 



(fc k 2 _)(kQ k 2 _) 



and the gap equations are 



1 
1 

AG 2 



[d 2 k] 
[d 2 k] 



fmiiko, fci;mi) 
(fc 2 , — k 2 _){k^ — fc 2 ) 

f M 2(k ,k 1 ;m 1 ) 
(fcg — k 2 _){k,Q — fc 2 ) 



(11) 
(12) 

(13) 
(14) 



where [d d k] = d d k/{2^) d , and *4(mi) = a{m 1 )±2b{m 1 ), 
with 

a(m) = k\ + fj 2 + M 2 + m 2 , 

b(m)= [ M 2 fc 2 + (pi 2 + M 2 )m 2 ] 1/2 . 

We also have 

fml ( k o, fa ; m) = k 2 - k\ + (y 2 + M 2 - m 2 ) , 
f M 2(k Q , ki;m) = fcg - k\ - (p? + M 2 - m 2 ) . 

Next, we perform the fco integral. Generically, we have 

f(ko,fa;m) 



[dfa 



[dfa] 



[dk ] 

, 1 
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(fc k 2 _)(k'^ fc 2 ) 
/(fc+,fci;m) f(k^,ki;m) 



(15) 



fc + 
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Hence, the gap equations become 



temperature effective potential corresponds to 



2gi 



AG 2 



-A 



r l 
K 
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~kZ 
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(fc^ fel)]' 7«(mi ) 0,M) = -y A [dfe 1 ]{fc+(mi) + fc_(mi) (18) 



m 

~6 Vfc 



M 2 / J 
1 



+ 



The fci integrals are logarithmically divergent and need to 
be regularized by imposing a cutoff A. The gap equations 
can be obtained by direct differentiation of the effective 
potential, V e s(mi,0, M), corresponding to 

r A r i 

V (1) (mi,0,M) = -/ [dfci] fc+Oi) + k-(m 1 ) . (16) 

The finite temperature gap equations can be obtained 
formally from the zero temperature ones, with the re- 
placements 



±^±.[l-2n F (k-)], 



(17) 



where np(k) = [e /3fc + 1] 1 is the Fermi-Dirac distribu- 
tion function, with j3 = 1/T. It follows that the finite 



| Mi 



(mx)] +-ln[l 



K)]} 



B. asymmetric case: rri2 7^ 

We have [H: h T h = A + B ■ a + (C + D ■ a)r 3 and 
/i/i T = A' + B' •<?+ (C + 3' ■ ct)t 3 . It is convenient to 
use an explicit isospin matrix representation, i.e. 



h T h 
hh T 



A + C+(B + D)-cr 

A-C+{B-D)-a 

A'+ C'+ (B'+ D') ■ a 

A'-C'+(B'-D')-(j 



By inspection, in the asymmetric case we find that the 
isospin trace gives two contributions similar to the sym- 
metric case, corresponding to masses m± — mi ± m^- 
We have 



1 f A r 

V {1) = — / [dfciK k+{m + ) + fc_(m + ) + fc + (m_) + fc_(m_) 

2 Jo 



(19) 



\ ln[l + e -^+( m +)l +| ln[l + e-^-( m +)l +4 ln[l + e-^+( m -)l +| ln[l + e-^-( m -)l \ 
p a a a > 



2 





III. RENORMALIZED EFFECTIVE POTENTIAL 

It is sufficient to perform renormalization at zero temperature and chemical potential (T = // = 0). In this context, 
we can define the renormalized coupling constant in terms of the physical scattering amplitude of fermions at a 
particular momentum scale. The addition of fi and T will only result in finite corrections to the gap equations, and 
therefore to the vacuum values of m, and M . 



A. T = M = 



In order to expose the divergences in the k\ integral, we consider now the integral 



dk\/k 2 + m 2 = -\Jk 2 + m 2 + — \\\{k+\/k 2 + m 2 ) 



which indicates that we have both logarithmic and quadratic ultraviolet divergences. First, we eliminate the quadratic 
divergence by adding a constant term at a fixed momentum cutoff, A, i.e. 



V {1) {m 1 ,m 2 ,M)- 



1 



[dfci] fc + (m_) + fe_(m_) + k + {m + ) + fc_(m+) — 4fcj 



(20) 
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One more integral |2lJ yields the unrenormalizcd effective potential at T = \x = 0, as 

1 



V cS {m u M) 



= M2 {-h 



4tt 



(21) 



1 

8^ 



x9 2A , x9 2A x9 2A , x9 2A 

(M + ?ti+) 2 In — + (M - m+) 2 In — + (M + m_ ) 2 In — + (M - m_) 2 In - 



|M + m+| 



\M-m+\ 



\M + m-\ 



\M-m-\ 



Noting that 2(M 2 + m\ + m\) = (M 2 + to , ) + (AI 2 + m 2 _ ) , we renormalize by requiring that the renormalized coupling 
constants, g 2 and G 2 , satisfy 



Mo 



9i.R 



mMo CP- 



(22) 



Here, the masses {niiQ,Mo} have arbitrary renormalization values on which the coupling constants will depend. 
Eqs. |2"2"|) are solved for g 2 and G 2 as a function of g 2 R and G\ . We obtain 



and 



1 

4G 2 



1 

2^F 



47 ~4^ n (2A) 5 



2<4 



_L J_i _L_ 

2tt 4tt n (2A) 2 



1 m 



16tt Mo 



hi 



M - m+o 



M + m+o 



1 m_ 



16tt Mo 



In 



M - m_ 



Mo + m_o 



where to±q = m w i m2o , and we have introduced the renormalization scale 



7 2 =(M 2 - ? ti 2 - m 2 ) 2 - 4m 2 m 2 = (AI 2 - m%)(M 2 - ? ti 2 _ ) 



(23) 



(24) 



(25) 



Hence, the renormalized effective potential at T = fj, = is obtained by substituting in Eq. I)21|l. the feare coupling 
terms as given by Eqs. I|23|l and (|24|l . We obtain 



Kff(TO*,M) 



M 2 /3 + m\oi\ + m\a.2 + 
Mm+ 



Ar- 



ln 



| M 2 - m\ | 



A/ 2 



8tt 



In- 



4tt 



In 



M - m+ 


Mm_ 


In 


M - m_ 


M + m+ 


47T 


M + TO_ 



(26) 



with 



2<& 



3 

47T 
1 

2^ 



1 m 



16tt M 



In 



M - m+o 



M + m+o 



1 m_ 



16tt M 



In 



M - m_ 



M + m_ 



Correspondingly, the gap equations are 



a x 1 , M 2 

— H In 

2 8tt 



7 



on 1 , Af 2 - m 2 

— H In- 1 

2 8tt 



n 2 



m + ( 
+ V 2 



1 |M 2 - 
m 

8n 7 



In 



|M 2 -m 2 | 1 



7 



8tt 



In- 



Q!2 

Y 

M 2 - ?n 2 

7 



87T 7 



87rV 1 



AI - 



77l-| 
"8^ 



hi 



M 
"8^V 
M - m + 



fin 



M + 771+ 
M — 771+ 



AI + m + 



+ln 



-In 



AI - m_ 



M + 771- 
M - 771- 



AI + m. 



M + tbj 



8tt 



ln 



M - 77l_ 



M + m_ 



= 0, 



with 



1 

47T 



£ = + 



4tt 



= 
= 



(27) 
(28) 

(29) 
(30) 
(31) 

(32) 



The solutions to,; = m* and M = M* give the local extrema of V e e, and represent physical parameters that must 
be independent of the renormalization scale 7. We note that if we solve for the combinations Si = j3 — oti, the 
renormalization scale 7 drops out. Therefore, Si are true physical parameters in the theory, and their values control 
which of the condensates and M can exist. The third physical parameter in this model is AS = S% — 81 = a\ — 01%, 
which is also independent of 7. Of course, only two of the three physical parameters Si, S2 and AS, are independent 
of each other. 
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B. Finite T and fi 

As advertised, the subtractions necessary to remove the ultraviolet divergences at T = fj, = 0, are sufficient to 
renormalize the effective potential at finite temperature and chemical potential. In order to see this, we note that we 
can write the bare coupling constants (|23[) and l|24|) as 



1 



a t + X 



1 

4G 2 



P + X 



(33) 



where X is the divergent integral 



y/k\ + (M + m+o) 2 ^k 2 + (M - m+o) 2 ^ + Wo + ™-o) 2 \A 2 + (M - m_ ) 2 
1 (2A) 2 

= - — In h terms that vanish as A — > oo . 

47T 7 

Then, the full renormalized effective potential at finite temperature and chemical potential can be written as 

1 



V cS {mi,M) = M 2 (3+mfai +m 2 2 a 2 - 



[dfci ] <^ k+ (m+ ) + (m+) +fc + (m_ ) + k- (m_) — 4fci 



(34) 



(35) 



4 ln[l + e-^+( m +)l +| ln[l + e-' 3fe -( m +)l +| lnh+e-^" 1 -)] +^ lnh+e" 

p p p p 

7l/ 2 + m?+m| 



/3fc_(m_)l 



Vfc 2 + (Mo + m +0 ) 2 ^k 2 + (M - m +0 ) 2 v^i + + m_ ) 2 v^l" + ( M o - m_ ) 2 

I 



IV. VACUUM 

At T = /x = 0, the minimum effective potential takes 
the form 



Kff (m^ M) = - ^-(M 2 + m 2 + m\) 

47T 

= - — (M 2 + m 2 +) - — (M 



(36) 



8tt 



8tt 



2 i 2 

1 m 



where, for simplicity, we have dropped the * notation 
of the parameter values at the minimum. We analyze 
now the solutions of the gap equations and determine 
the true vacuum of the theory as the global minimum of 
the effective potential. The parameter space corresponds 
to the triangle depicted in Fig. ^ Here, the corners of 
the triangle correspond to the situation when two masses 
are zero: 



• mi 



m 2 = 0: M 2 = 7e -(i+^/3) 5 an d 



_ _J_ -(1+4*0) 

4tt 

= 7e-( 1+4mi >, and 
F eff (m l! 0,0) = --I- 

47T 

M = 0: m\ = 7e -(i+47ra 2 ) ! and 
V/ off (0,m 2 ,0) = -^ e -(i+*™»). 

47T 



• mi 



Kff(0,0,M) 



• M = m 2 = 0: m\ 



(37) 



(38) 



(39) 



m o =0 




FIG. 1: (Color online) Parameter space for the effective po- 
tential at T — fi = 0. The corner with M — m,2 = indicates 
the QCDdike sector. 



while the sides of the triangle correspond to the case 
where one mass is zero. We have 

• M = 0. Here, we have m+ ^ m_ (m 2 ^ 0), and 
the gap equations become 



( - 1 i m 
mi ct\ H m — 

V 47T 



m 2 

47T 



In 



mi — m 2 



1 |m 
m 2 ( a 2 + — m — 

47T 



mi 
In 

47T 



mi + m 2 



mi — m 2 



mi + m 2 



= 0, 
(40) 

= 0. 
(41) 
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TABLE I: Position of the global minimum of the effective 
potential at T — fi = 0. 



Parameters 



masses 



V c s(mj,M) 



Si < 8 2 < (AS > 0) 
8 2 < Si < (AS < 0) 
Si < < 8 2 (AS > 0) 
< Si < 8 2 (AS > 0) 
< 82 < Si (AS < 0) 
6 2 < < <Ji (AS < 0) 



M > mi > 7712 
M > 7712 > mi 
7772 > M > 777i 
7772 > m-i > M 
777i > 7772 > M 
7771 > M > 7772 



V«ff(0,0,M) 

K eff (0,0,M) 

Kff(0,m 2 ,0) 
V eff (0,ro 2 ,0) 

v; ff (mi,o,o) 

Kfi(mi,0,0) 



• mi = 0. In this case, wc have ?7i + = — m_ 
The gap equations lead to 



m 2 \a 2 



47T 7 



itr, 



M 

47T 



■In 



M 



m 2 



\ Air 



m 2 , 
In 



M + m 2 



M — m 2 



M + m 2 



= m 2 . 

= 0, 
(42) 

= 0. 
(43) 



When only one mass (mi, m 2 , or M) is zero, (see ap- 
pendix of Rcf . Q) the "global" minimum, along the side 
of the triangle, corresponds to one of its ends and which 
end exactly has the lower potential depends on the value 
of the parameters Si, d 2 , or AS, respectively: 

• If Si < 0, then V; ff (0,0,M) is the "global" mini- 
mum along the m 2 = line, whereas if Si > then 
the "global" minimum of the effective potential is 
V e g (mi, 0, 0). The critical value <5i = corresponds 

tO M = 777l. 

• If S 2 < 0, then V; ff (0,0,M) is the "global" mini- 
mum along the mi = line, whereas if S 2 > then 
the "global" minimum of the effective potential is 
V e s(0, m 2 ,0). The critical value S 2 = corresponds 
to M = rn 2 . 

• If AS = OLi — a 2 = S 2 — Si < 0, then V e e(mi,0, 0) 
is the "global" minimum along the M = line , 
whereas if A<5 > then the "global" minimum of 
the effective potential is V o s(0,m 2 ,0). The critical 
value A(5 = corresponds to mi = m 2 . 



• m 2 = and we recover the symmetric-case results 
with (mi ^ 0, M ^ 0): 



1 \M 2 -m\\\ M 



mi I cti + — In I — —In 

V 47r 7 



M(p + — In 

V 47T 



1 \M 2 - m\ 
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47T 



mi 

47T 



In 



M - mi 



M + mi 



M — mi 



M + mi 



(44) 



0. 



(45) 



Following an approach similar to the one described in 
the appendix of Ref . Q , one can show that the effective 
potential in the case when all three masses mi and M 
are nonzero, has a local minimum intermediate between 
the "corner" values. In conclusion, the global minimum 
of the effective potential at T = /j = is always located 
at one of the corners of triangle depicted in Fig. ^ bi 
Table [I] wc show which corner corresponds to the global 
minimum of the effective, as a function of the relative 
values of the parameters Si and 5 2 . 
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PHASE STRUCTURE 



We discuss now the effective potential ()35(l corresponding to the three possible vacuum phases identified in the 
previous section. We have: 

• mi = m 2 = 0. In this case our model reduces to the pure Cooper-pairing model with M* 2 = A 2 is the 
dynamically generated gap. Here, we choose mm = m.20 = 0, Mq = A 2 = 7, and f3 = 0. Then, we can write 

Veff(rrii,M) = m 2 ai +m 2 a 2 



(46) 



[dki]\ k + (m + ) + fc_(m + ) + fc + (m_) + fc_(m_) - Aki - 



2(M 2 +ml+mj) 
^Jk\ + A 2 



+ - D ln[l + e-^+( m +)l +| ln[l + e -^-(™+)l +| ln[l + e-' 3fe +( m -)l +| ln[l + e^ fc -( m -)] 
p p p p 

If we set mi = m 2 = 0, then the chemical potential becomes irrelevant and can be transformed away. We obtain 

2 



V oS (m t ,M) 



M 2 r M 2 
4tt L A 2 



1 



P Jo 



[dki] ln(l 



-0E\ 



lnfl + e"^) 



(47) 



with E = y/k\ + m\. 

The effective potential for the Cooper pair sector at finite temperature is fi independent. The model has a 

second-order phase transition to the unbroken phase at a critical temperature T c = (A/7r)e 7s , where 
je = 0.577 ■ • • is Eulcr's constant. 
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• M = m,2 = 0. This is the Gross-Neveu sector |24|. Here, m\ = mp is the dynamically generated fcrmion mass 



We choose M = m 2 o 



no 

,2 - 



7. Furthermore, we have g 2 R = t and «i = 0. Thus, we can write 



V e s(rrii, M) = M 2 (5 + m 2 2 a 2 
1 



(48) 



[dki]< k + (m + ) + fc_(m + ) + fc + (m_) + fc_(m_) - 4fci - 



2(M 2 



ml) 



+ | lnri + e-^+( m +)l + | ln[l + e-^-( m +>l + | ln[l + e -' 9fe + (m - ) l +| \n\l + e~ 0k - {r 
p p p p 



If we set M = m 2 = 0, then we obtain m+ = m_ = mi, fc± = E ± [i with = y^i + m i) an d the effective 
potential in the Gross-Neveu sector can be written as 






2 

r mf 






In —i - 1 




47T 


L mf 


/?7o 



> [dfc 1 ]{ln[l+e- /3 (' B +' 1 )] + ln[l+e-^ E -ri] } 



(49) 



We note that the effective potential in the Gross-Neveu sector is similar to the effective potential in the Cooper 
pair sector, with mp replacing A, but depends on the chemical potential \i at T =/= 0. 

The Gross Neveu model has spontaneous symmetry breaking at zero chemical potential and zero temperature. 
At zero temperature, the model undergoes a first-order phase transition to the unbroken symmetry phase as 
we increase the chemical potential [3> H3- At zero chemical potential, the symmetry is restored through a 
second order phase transition at the critical temperature @, T c = (mp/n) exp[-jE — 7/i 2 C(3)/(47r 2 T c 2 )] . Thus, 
the phase diagram of the Gross-Neveu model has a tricritical point, with approximate values |24j fi c /m F = 
0.608, T c /m F = 0.318. 

• mi = M = 0. Here, m* 2 2 = 5m 2 is the dynamically generated mass asymmetry. This case is identical with 



the Gross-Neveu sector case, with mp being replaced by 8m: We choose mio = M 
9ir = n ( or "^2 = 0). The effective potential becomes 



0, m 



2(1 



V eS (mi,M) = M 2 (3 + mfaj 



(5m 2 = 7, and 
(50) 



[dfciK k + (m + ) + fc_(m + ) + fc + (m_) + fc_(m_) - 4fci 



2(M 2 



ml) 



I ln[l + e-^+( m +)] +\ ln[l + e-" fc -( m +>] +- ln[l + e- 3fc +( m -)] +- ln[l+e 



y/kY+Sm? 



P 1 P 

When mi = M = 0, we obtain 
V cS (m u M) 



ml 


2 

, mi 






In— | - 1 


--/ 


47T 


- m F 





[dfci]{ln[ 



1+e 



ln[l- 



-/3(S-P) 



]}• 



(51) 



with £" = Y/fcf 



VI. CONCLUSIONS 

In this paper we studied a generalization of a simple 
model with two spontaneously broken symmetry phases 
first discussed in Ref. The model is studied in (1+1) 
dimensions within the leading order in large N approxi- 
mation. We find that the phase diagram with two differ- 
ent masses is similar to the case when the masses are the 
same in that depending on the values of the two renor- 
malized coupling constants one is always in one of three 
distinct phases. 

A generalization of this model to higher dimensions in- 
cluding the case of spatially inhomogeneous order param- 



eters is under development. We expect that the results 
found here in 1+1 for homogeneous order parameters will 
persist in mean field regardless of the number of dimen- 
sions. Thus, in order to obtain phase coexistence we will 
have to consider inhomogeneous condensates. For this 
type of condensates we expect a small coexistence region 
similar to the Larkin-Ovchinikov-Fulde-Ferrel state [8( in 
ferromagnetic superconductors, while preserving the tri- 
critical point. 

After this work was completed, we became aware of a 
new study of the massive Gross-Neveu model using semi- 
classical methods in the large-N limit [2(J. These semi- 
classical methods display a kink-antikink crystal phase 
and it would be interesting to also perform such an anal- 
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ysis for the model of this paper. In 1+1 dimensions, semi- 
classical analyses such as those described by Dashen, 
Hasslacher and Neveu [2^ often gave exact answers for 



the S-matrix (especially in exactly solvable models) and 
include more information than the straight forward study 
of the effective potential in the largc-N approximation. 
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